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Abstract 

A full (3+l)-dimensional calculation using the Lagrangian hydrodynamics is proposed 
for relativistic nuclear collisions. The calculation enables us to evaluate anisotropic flow 
of hadronic matter which appears in non-central and/or asymmetrical relativistic nuclear 
collisions. Applying hydrodynamical calculations to the deformed uranium collisions at 
AGS energy region, we discuss the nature of space-time structure and particle distributions 
in detail. 
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1 Introduction 

The study of the hot and dense matter which is produced in relativistic nuclear collisions 
has received the intensive attention, and a number of experiments have been done to 
realize the phenomena ||]. Hydrodynamical model is one of the established models for 
describing global features in the relativistic nuclear collisions as collective flow. Recently, 
anisotropic flow phenomena have been observed at AGS [§, [j| and also at SPS ||, ||, |(|. 
Several authors have argued the relation between the behavior of the collective flow and 
the equation of states. Based on the relativistic hydrodynamical model, Rischke reported 
that the existence of the minimum point in the excitation function of directed flow would 
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suggest the phase transition |7j. Danielewicz has shown that the elliptic flow is sensitive 
to the difference of the nuclear equation of state by using a relativistic hadron transport 
model [pi. Sorge has discussed centrality dependence of elliptic flow based on the event 
generator RQMD which includes phase transition 0, [|. Hence, analysis of the flow can 
inform us the property of the nuclear equation of state and quark-gluon plasma (QGP). 

Since Bjorken proposed the scaling solution 1C], a number of investigations based on 
the relativistic hydro dynamical model have been done. Those reproduced successfully 
experimental inclusive spectra at both AGS and SPS energies 1{17[|. Ornik et al. in- 
vestigated single particle distribution and Bose-Einstein correlation at SPS Sollfrank 
et al. discussed the hadron spectra and electromagnetic spectra at SPS IE]. Hung and 
Shuryak mentioned the equation of states, radial flow, and freeze-out at AGS and SPS 
JT(|. Morita et al. discussed single particle distribution and Bose-Einstein correlation at 
SPS [ 18 1 . Therefore the relativistic hydro dynamical model works well in the analyses 
of many kinds of phenomena in the ultra-relativistic collision. However, in most of the 
studies based on the relativistic hydrodynamical model, cylindrical symmetry is assumed 
and, therefore, discussions are limited only to the central collisions. 

Recently anisotropic flow has been analyzed by using relativistic hydrodynamical 
model. Teaney and Shuryak proposed "nutcracker" phenomenon and Kolb et al. 
discussed anisotropic flow and phase transition | |20| |. But these analyses use Bjorken's 
scaling solution in the longitudinal direction, and their discussion is restricted in the 
mid-rapidity region. 

In order to investigate anisotropic flow, the quantitatively reliable (3+l)-dimensional 
relativistic hydrodynamical calculation is indispensable. The full (3+l)-dimensional cal- 
culation of relativistic hydrodynamical equation has already been done by Rischke et al. 
1 2 1 ] and Brachmann et al. [22]. Rischke pointed out that the minimum of the excitation 
function of directed flow suggests the existence of the phase transition |2l|] . Brachmann et 
al. [p2} | discussed the antiflow of nucleons at the softest point of the equation of state using 
three-fluid dynamics. Their numerical schemes are based on the Eulerian hydrodynamics. 

Here, we present the Lagrangian hydrodynamic simulation without assuming cylindri- 
cal symmetry which makes the full (3+l)-dimensional analyses possible. The Lagrangian 
hydrodynamics has several advantages over the Eulerian hydrodynamics to investigate 
phenomenon in ultra-relativistic nuclear collisions. As is well known, in Lagrangian 
method, grid points move along the flow, and is superior to Eulerian when calculation 
region varies rapidly. This is a great advantage for high energy collision where initial 
nuclei are Lorentz-contracted. Secondly, we can follow a trajectory directly in the phase 
diagram in case of Lagrangian algorithm. This allows us to study heavy ion collision 
phenomena together with the equation of state. 

As an example of the anisotropic flow in heavy ion collisions, we apply our hydro- 
dynamical model to deformed uranium-uranium collisions in AGS energy and analyze 
the elliptic flow in detail. Shuryak pointed out the remarkable features of the deformed 
uranium collisions which are suitable for the important problems such as hard process, 
elliptic flow, and J/tp suppression [23|. The effect of such deformation to the flow [24] and 
J suppression p5| has been investigated. Here, we discuss the influence of deformation 
of uranium nucleus on elliptic flow. 

In Sec. II we introduce relativistic hydrodynamical equation and explain our original 
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algorithm of the numerical calculation of (3+l)-dimensional relativistic hydrodynamical 
equation. In Sec. Ill we apply our hydrodynamical model to investigation of elliptic flow 
which is produced in deformed uranium collisions and discuss the effect of the deformation 
on elliptic flow in detail. Section IV is devoted to the summary of the paper. 

2 Algorithm to Solve Hydrodynamical Equation 
2.1 Relativistic Hydrodynamical Equation 

The relativistic hydrodynamical equation in Lorentz-covariant form is 

d^ v = 0. (1) 

Since we discuss a system formed by nuclear collision, baryon number current conservation 
should be also taken into account, 

dull = 0- (2) 
In this paper, T^ v is taken as energy-momentum tensor of the perfect fluid, 

= eu»u v - p{jf - u»u u ), (3) 

and baryon number current is given by 

f B = n B {T,ii)u v . (4) 

Here e, p and n B are energy density, pressure and baryon number density, respectively. 
These are the functions of the coordinates through the temperature T{xu) and the baryon 
number chemical potential = 7(1, v x , v y , v z ) and g^ v = diag.(l, —1, —1, —1) are 

local four-velocity and metric tensor, respectively. If the equation of state is properly 
given, we can solve the coupled eqs. (p]) and (^) and obtain the chronological evolution 
of temperature and chemical potential. In order to make our numerical method clear, 
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eqs. (Q) and (Q) are rewritten as 
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where 7 = l/VT" 

eq. one can obtain entropy conservation law, 

3% = 0, 

with the aid of the thermodynamical relation, 

e +p = Ts + /m B , 



(5) 

e + p. From time-like projection of eq. (|l|), u^d^T^ = 0, and 



(6) 
(7) 



where s^=su^ is the entropy current density. We solve numerically eqs. @ and (||) with 
the algorithm which will be explained in the next subsection. 



2.2 Computational Scheme 

Most hydrodynamic calculations which are used for investigating various phenomena in 
heavy ion collisions are based on the Eulerian hydrodynamics. Sollfrank et al. analyze 
hadron and electromagnetic spectra by using SHASTA algorithm . HYLANDER and 
HYLANDER-C algorithm are used by Ornik et al. jTJ] and Schlei and Strottman Q, 
respectively. Rischke et al. use RHHLE algorithm and study hydrodynamics and collective 
flow §. 
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Here, we solve the (3+l)-dimensional relativistic hydrodynamical equation with the 
Lagrangian hydrodynamics. The Lagrangian hydrodynamics has several advantages over 
the Eulerian hydrodynamics to treat ultra-relativistic nuclear collisions. At high energies, 
initial distribution of energy localizes due to collision of Lorentz contracted projectile and 
target. To treat the situation, fine resolution is required in the Eulerian hydrodynamics 
and computational cost becomes expensive. On the other hand, in the Lagrangian hy- 
drodynamics, discretized grids move along the expansion of the fluid, therefore we can 
perform the calculation at all stages on the lattice points which we prepare in the ini- 
tial condition. For example, in our previous calculation the fluid expands four times 
larger in longitudinal direction. This fact means we need four times number of grid if we 
use naive Eulerian type algorithm. Another merit of the Lagrangian hydrodynamics is 
that it enables us to derive the physical information directly, because it follows the flux of 
the current. For example, the path of a volume element of fluid in the T-fi plane can be 
traced as we will demonstrate in the next subsection. Therefore, we are able to discuss 
how the phase between hadron phase and QGP phase effects the physical phenomena by 
the Lagrangian hydrodynamics. 

Our numerical calculation of the (3+l)-dimensional relativistic hydrodynamical equa- 
tion is as follows: First, the coordinate at time t + At, x m = X m (t,i,j,k) is evolved 
as, 

X m (t + At, i,j, k) = X m (t, i,j, k) + ^jhhll^- At. (8) 

U [T, l, J, K) 

By definition of the Lagrangian hydrodynamics, the coordinates move in parallel with j M 
and s^. 

Second, the local velocity is determined, 

v m {t + At,i,j,k) = v m (t,i,j,k) +d t v t {i,j,k,t)At 

3 

+ J2d n v m (i,j,k,t)(X n (t + At,i,j,k) - X n (t,i,j,k)), (9) 

n=l 

where the value of d^v^ is obtained from eq. (||). 

Finally, the temperature and chemical potential are derived. The volume element da 11 
at time t is surrounded by eight points, X^it, i, j, k), X^{t,i + l,j,k), X^(t,i,j + l,k), 
X^(t,i,j,k + 1), X' J -(t,i+ l,j + 1, k), X^{t,i + 1, j + 1, k + 1). Using this volume 
element, eqs. (|2|) and @ are rewritten as 

s(T(t + At, i, j, k),fx(t + At, i,j, k))u\t + At, i,j, k)da t (t + At, i, j, k) 

= s(r(t,i,j,fc),/x(t,i,j,fc))u i (t,i,j,fc)do-j(t,i,j,A;), (10) 

n B (T(t + At , i , j, k) , n(t + At , i , j, k) )u* (t + At , i, j, k)da t (t + At,i, j, k) 

= n B (T(t,i,j,k),n(t,i,j, k))u l {t,i, j,k)da t (t,i, j,k). (11) 

Here, by virtue of the determination of coordinates eq. (|8|), we can use the relation, 

u^da^ = u dat. 
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Since s and n B depend on T and fx, using up to the first order of the differences of 
temperature, AT(t, i,j, k) = T(t + At, k) — T(t, k), and of the chemical potential, 
A/i(t, i, j, k) = fi(t + At, i, j, k) — fi(t, i,j, k), we expand s and n B as, 



s(T(t + At, i, j, k),n(t + At, i, j, k)) 
= s{T(t,i,j,k),fi(t,i,j,k)) 
ds 
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Substituting eqs. (12) and ( |i~3| ) to eqs. (10) and (0), we obtain the temperature and 



chemical potential at the next time step, 
T(t + At,i,j,k) = T(t,i,j,k) 
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where A Sjns and 
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respectively. These numerical procedures are the extension of the method in ref. [|l^] . 

In this algorithm CPU time is almost proportional to the number of the lattice point. 
Numerical calculation of the relativistic hydro dynamical equation in this paper has been 
performed at the Institute for Nonlinear Sciences and Applied Mathematics, Hiroshima 
University. Average floating point operations for (55, 63, 49) space points and 5300 time 
steps are 36 Tera in each calculation reported in the next section. 

In order to analyze anisotropic flow with high accuracy, the artificial anisotropy which 
can be caused by the discretization of space should be small. We checked the reliability 
of our calculation by comparing the results with rotated spatial grid. Figure |] shows that 
results of flow obtained by different choice of the grid. The difference of flow between two 
results is less than 0.15 % in the present calculations. 
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Figure 1: The left part of the figure indicates the flow at y ~ fm and t = 20 fm/c in Au + Au 
10.8 AGeV collision [p6|| . The right part of the figure shows the same result which is obtained 
by the calculation with the grids rotated 26 degree on the x-z plane. 



2.3 Path in the Phase Diagram 

The Lagrangian hydrodynamics enables us to trace easily the history of the trajectory 
of flux. In order to test the applicability of our algorithm to study the chronological 
trajectory of the volume element in the phase diagram, we use the equation of state 
which contains the first order phase transition only in this subsection. Above the phase 
transition, the thermodynamical quantities are assumed to be determined by QGP gas 
which is dominated by massless u, d, s quarks and gluons. In the QGP phase the pressure 
is given as, 



(32 + 21^ 4 iV//Mf r 2 , N f 
1 180 2 V3/ 4vr 2 



B, 



(16) 



where Nf is 3 and B is Bag constant [15, 16]. For the hadron phase we use the excluded 
volume model [27| which contains all resonances up to 2.0 GeV |28|| . In the hadron phase 
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the pressure for fermion is given as, 



p exd (T,{^}) = ^pf^iT^i-Vap^iT,^})), 

i 



(17) 



where p ldeal is the pressure of ideal hadron gas and Vq is excluded volume of which radius 
is fixed to 0.7 fm. Putting the critical temperature as 160 MeV for zero chemical potential, 
the Bag constant, B*, is given as 233 MeV. Figure [2] shows the equation of states as a 
function of temperature and chemical potential. Figure [3| indicates the phase boundary 
which is determined by the pressure balance between the two phases, i.e., pq = pn- In the 
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Figure 3: The phase boundary which 
is determined by the pressure balance. 



mixed phase we introduce the fraction of the volume of the QGP phase, \{xn) (0 < A < 1) 
and parameterize energy density and baryon number density as, 

e M (X,T*(fj,)) = Ae Q (T*( M ))-(l-A)eff(T*( M )), 
n B M(X,T*(fj,)) = XnBQ(T*(pL))-(l-X)n B H(T*(fi)), (18) 

where T*((i) is the value of temperature on the phase boundary in phase diagram. Con- 
trary to in the Eulerian hydrodynamics where the boundary condition should be consid- 
ered on the discontinuous plane between two phases |2^, |3(| , in our algorithm, by virtue 
of the explicit use of the current conservation equations, the flux of the fluid can be traced 
easily even if the discontinuity of the first order phase transition exists |l8| , 31]. For initial 



conditions, we use the results which are obtained by URASiMA (Ultra-Relativistic A- A 
collision Simulator based on Multiple scattering Algorithm) in Au+Au 20 AGeV collisions 



[26]. 



Figure |] shows the typical paths in the phase diagram. For instance, the trajectory 
of the volume element of the grid number (x,y,z) = (0,-1,0), which starts from QGP 
region, moves along phase boundary before entering the hadron phase. The volume ele- 
ment of (0, —5,0), which starts from the mixed phase, also moves along phase boundary 
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Figure 4: The paths in the phase diagram. The dot-dashed line stands for the path which 
starts from the QGP phase and on this line Ub/s is 0.065. The solid line stands for the path 
which starts from the mixed phase, and on this line %/s is 0.059. The dotted line stands for 
the path which starts from the hadron phase, and on this line ns/s is 0.044. 



before turning into the hadron phase. On the other hand the volume element of the grid 
number (x,y,z) = (0,-7,0), which starts from hadron phase, draws smooth trajectory 
to the freeze-out. This behavior is the same result as discussed in ref. [16|. By virtue of 
the Lagrangian hydrodynamics we can easily trace the trajectory which corresponds to 
the adiabatic paths in the T-\x plane. 



3 Effect of Deformed Uranium on Collective Flow 

As an application of our hydro dynamical model, we investigate the flow which is generated 
in the deformed uranium collisions. Shuryak has pointed out that remarkable problems 
such as hard processes, elliptic flow, and the mechanism of J ftp suppression can be resolved 
by using the deformed uranium collisions [23]. Since Danielewicz showed that elliptic flow 
is sensitive on nuclear equation of state in AGS energy ||], the elliptic flow is one of the 
hottest topics in high energy nuclear physics. The high accuracy experiments for elliptic 
flow have been done at AGS |2|, fj] and SPS ||, |6| . Recently, using the cascade model 
of the ART, Li discussed the elliptic flow in the deformed uranium collisions. If the 
deformation causes a large influence on anisotropic flow, the analyses of collective flow 
using U+U collisions are promising for investigating the difference between QGP states 
and hadron states. In order to analyze anisotropic flow of deformed uranium collisions, the 
(3+l)-dimensional relativistic hydro dynamical model plays a central role, i.e. it provides 
us reliable quantitative results. 

Here, the ellipticity is measured by the asymmetry of azimuthal distribution of particle 
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which is expanded based on Fourier series, 



dN 

— ~ vq(1 + 2vi cos{4>) + 2v 2 cos(20)), (19) 
a<p 

where <ft is azimuth and Vq is normalization. Parameters v± and v% correspond to the 
intensity of the directed flow and elliptic flow, respectively. 

The shape of deformed uranium nucleus is approximately ellipsoid and short (Rt) and 
long (R c ) semi-axis are given as, 

Rt = R s (l-^S), 

Rc = R S (1 + ^S), (20) 



where <5 =0.27 is the deformation parameter [32]. We will investigate how the deformation 



and the orientation between two colliding deformed uranium nuclei influence the flow of 
produced hadrons. Among many types of collisions for the orientation, we focus on two 
types of collision, i.e., tip-tip collision in which the long axes of two nuclei are along beam 
direction, and body-body collision in which the long axes of two nuclei are parallel each 
other but perpendicular to the beam direction. We calculate also sphere-sphere collision 
for comparison. 

3.1 Model Description 

The outline of our calculation procedure is as follows: First, we parameterize the initial 
conditions of energy density, baryon number density and local velocity based on the result 
of event generator URASiMA. 



Our event generator URASiMA p3|, 34, pq] is characterized by multi-chain model 



(MCM) by which multi-particle production process can be described successfully. In 
URASiMA the detailed balance between quasi-two-body production and absorption pro- 
cesses holds. It is applicable to AGS and SPS energy regions and calculated results 
reproduce experimental data of hadron spectra |3^] . Recently thermodynamical proper- 
ties of hot and dense hadronic gas are also investigated by URASiMA pL E$5|. For 

a more 



detailed discussions, see refs. J33 



In order to solve the relativistic hydrodynamical equation, we need to introduce an 
equation of state. Since our calculation does not rely on any artificial assumption, we can 
investigate how the difference of equation of state has an effect on physical phenomena. For 
the first trial, we adopt the equation of state of the ideal hadron gas including resonances; 
this is the same equation of state as we used in the initial conditions. The temperature 
and chemical potential of volume elements vary with the space-time evolution of fluid 
until hadronization process occurs. 

Finally hadron spectra are obtained by Cooper- Frye formula |3(|. We assume that the 
hadronization process occurs when the temperature and chemical potential of the volume 
elements cross the boundary (solid line in fig. ||). The solid line is obtained so that 
the freeze-out temperature becomes 140 MeV at vanishing chemical potential, based on 
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Figure 5: The freeze-out condition. The solid 
line stands for freeze-out condition which is de- 
termined from thermal freeze-out (the dotted 
line) and chemical freeze-out (the dashed line). 



chemical freeze-out model and thermal freeze-out model ||37f |. Using Cooper-Frye formula 
the particle distribution is given as, 

E ^p = E ^3 / da ^eM(P^-^ f )/T f ]±V (21) 

where is degeneracy of hadrons and Tf and fit are the freeze-out temperature and 
chemical potential shown in fig. ||, respectively and is local velocity of the fluid on the 
hypersurface da^. We determine hypersurface da^ by evaluating the normal vector to the 
freeze-out hypersurface ||38| . We evaluate intensity of the elliptic flow V2 using eqs. ([l9|) 
and @. 

3.2 Initial Conditions 

In high energy collisions such as SPS and RHIC, the contribution of the spectator to ini- 



tial conditions is small and "in-plane" elliptic flow is enhanced [39|. However, if we focus 
on AGS energy regions in which incident energy is not so large, the effect of deformation 
can remain strong. In order to prepare appropriate initial conditions, we estimate the 
energy density distribution and baryon number density distribution by the event gener- 
ator URASiMA. We calculate the space-time evolution in U+U 10 AGeV and 20 AGeV 
tip-tip, body-body and sphere-sphere collisions at b = and 6 fm, using the event gener- 
ator URASiMA, where b is the impact parameter. In the case of tip-tip and body-body 
collisions, we consider the collision of the deformed uranium according to eq. (|2~0|). We 
assume that the hydro dynamical expansion starts when the projectile nucleus finishes 
passing through the target nucleus. At this time initial conditions of hydro dynamical 
model, T{Xfj), /^(x^), and u^(x^) should be given. The results of URASiMA are listed in 
table |(a) (U+U 10 AGeV) and table |(b) (U+U 20 AGeV). In body-body collision the 
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Table 1: The energy and baryon number densities obtained by URASiMA at (x, y, z) = (0, 0, 0) 
for different type of collision and the impact parameter, (a) U+U 10 AGeV, (b) U+U 20 AGeV 
Table 1(a): The result of URASiMA (U+U 10.0 AGeV). 



b [fm] 


type of collision 


initial time [fm/c] 


e [GeV/fm 3 ] 
at (0,0,0) 


Ub [fm 3 ] 
at (0,0,0) 





tip-tip 


8.5 


1.975 


0.977 


6 


tip-tip 


8.5 


1.403 


0.718 





sphere-sphere 


7.0 


1.807 


0.929 


6 


sphere-sphere 


7.0 


1.611 


0.809 





body-body 


6.5 


1.681 


0.819 


6 


body-body 


6.5 


1.490 


0.728 



Table 1(b): The result of URASiMA (U+U 20.0 AGeV). 



b [fm] 



type of collision 



initial time [fm/c] 



e [GeV/fm 3 ] 
at (0,0,0) 



ub [fm 3 ] 
at (0,0,0) 



tip-tip 



6.0 



2.871 



1.124 



tip-tip 



6.0 



2.060 



0.825 



sphere-sphere 



5.0 



2.744 



1.064 



sphere-sphere 



5.0 



2.241 



0.869 



body-body 



4.5 



2.248 



0.933 



body-body 



4.5 



2.253 



0.860 
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energy density and baryon number density are the smallest among different types of col- 
lision in the same incident energy and impact parameter. For tip-tip collision a decrease 
in the energy density and baryon density to impact parameter is the largest among the 
three types of collision. Then the initial energy density distribution and baryon number 
density distribution are parametrized based on these data in tables [j](a) and [j](b). The 
initial energy density and baryon number density distributions are given as, 

e = e max B(x,y,z), 
nB = n Bmax 5(x,y,z), (22) 

where B(x,y,z) is the distribution function, which is determined so that the result of 
URASiMA is reproduced, and e max and n^max are the values of the result of URASiMA 
in central region. We interpolate and/or extrapolate the values of e max and n# max from 
the data of tables |l](a) and [l](b) as follows: The distribution function is given as, 



^1 O sx " "sy & "sz 



(x + x cs ) 2 (y + y cs y (z + z cs ) z 

+ai exp(- 



(7^ (7^ (T^ 

w sx sy sz 

+ exp(- (x ~ 2 Xcp)2 - (j/ ~f p)2 - {Z ~ Z 2 cp)2 )}, (23) 
®px ®vy ®p z 

where c\ is normalization and the parameter, which corresponds to the ratio of energy 
density of spectator to that of participant, a\ is fixed to be 0.7 for all cases. The ratio 
of baryon number density is also fixed by ai. In eq. fl2"3| ) x cp , y cp , and z cp are the center 
of participant, and x cs , y cs and z cs are the center of spectator. Their specific values are 
determined geometrically by the position of projectile nucleus and target nucleus. For the 
participant, 

•^cp Ucp Z<:p 0) 

b 

&px — a3{Rx — 



R x V x 4 ' 

a 4 — . (24) 

7 
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For the spectator, 



%cs 
Vcs 

%cs 

o R7 , 



R x (b > R 2 

^ %-y/2RJ> - b 2 {b<R a 



Osz = a&-fr~- (25) 

R x ,R y , and R z are the radius of projectile nucleus and target nucleus in the x, y, and z 
directions. In eq. (|25|) ti is the initial time when the hydrodynamical evolution starts. In 
eqs. (24) and (|25|), parameters 03 ~ qq are tuned so that the energy density and baryon 



number density distribution, eq. (|22j), reproduce the result of URASiMA. The parameters 
<23 ~ ae are fixed to 0.7, 2.5, 1.3, 0.7, respectively. Figure ^ shows the energy density 
distribution of initial conditions in U+U 20 AGeV sphere-sphere collision at b = 4.0 fm. 
From this figure we can see clearly the contribution from participant part and spectator 
part. For the initial velocity distribution we neglect the transverse flow because its value 
is small in the result of URASiMA. The initial flow distribution is given as, 

v x (x,y,z) = 0, 
v y (x,y,z) = 0, 

v z (x,y,z) = -B v (x,y,z). (26) 
Here the distribution function B v (x,y,z) is given by, 

B v (x,y,z) = L [eM J^^_^l_^l ) 

Co (7 (7 (7 

°z u sx u sy u sz 

(x + x cs ) 2 (y + y C s) 2 (z + z cs ) 2 



+ exp(- 



<J 2 (J 2 (T 2 

" sx sy sz 



"cp) (U Ucp) 2 {Z Zap)" 2 , . . 

+a 2 exp( „— ^ „— ^ 2 )> ( 27 ) 



®px ®vv ®pz 



where c 2 is normalization and a 2 is parameter and fixed to be 1.4. Figure ^ shows the 
initial velocity distribution in the same case as fig. ^. The flow in longitudinal direction 
is similar to Bjorken's scaling solution in the central region. 



3.3 Calculated Results 

Figures [8| and |^ indicate the expansion of energy density and baryon number density in 
the central region. There is a slight difference in the life time in each case. The initial 
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Figure 6: The initial condition of energy density distribution at y = fm in U+U 20 AGeV, 
sphere-sphere collision. In this case the impact parameter is 4.0 fm. 
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Figure 7: The initial condition of velocity distribution at y — fm in the same condition as 

fig. |. 
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Figure 8: The time evolution of energy den- 
sity of central region in U+U 20 AGeV with 
6 = fm. The solid line, the dot-dashed line, 
and the dashed line stand for tip-tip collision, 
body-body collision, and sphere-sphere colli- 
sion, respectively. 
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Figure 9: The time evolution of baryon num- 
ber density distribution in the same case of 
fig. H The solid line, the dot-dashed line, 
and the dashed line stand for tip-tip colli- 
sion, body-body collision, and sphere-sphere 
collision, respectively. 



energy density and baryon density of body-body collision are the smallest among all three 
types of collisions and the difference between tip-tip collision and sphere-sphere collision is 
small. At initial time the difference of energy density among all three types of collision is 
Ae ~ 10 ~ 20% and the difference of time is At ~ 10 ~ 20% at e = 1 GeV/fm 3 . However 
those kinds of difference among collision types do not appear at final time. 

Figure 10 shows that behavior of V2 of nucleon as a function of the impact parameter 
for each type of collision at 10 and 20 AGeV. The elliptic flow parameter, V2, increases 
with impact parameter and reaches a peak at b = 3,4 fm and decreases in every case. 
Furthermore for body-body collision vi does not vanish at b = 0. Consequently, the effect 
of deformation is not negligible. In order to make this characteristic behavior clear we 
focus on the pressure distribution. 

Figure 11 displays the pressure distribution in U+U 10 AGeV sphere-sphere collision. 
From the figures at z ~ fm we can see that the pressure gradient in the y direction 
increases with the impact parameter. Therefore one might consider that i>2 increases with 
impact parameter, because the velocity of produced particles in x direction is larger than 
one in y direction. But fig. 10 shows that the value of V2 starts to decrease at about 
b = 3 ~ 4 fm. Here we focus on the pressure distribution in z-plane in fig. 11, and we 
can see that the effect of spectator increases with the impact parameter. Because the 
spectators block the flow in x direction, the growth of flow in x direction is suppressed. 
Consequently the behavior of V2 is determined by both of the pressure gradient and the 
effect of spectator. 

Figure ^ shows the pressure distribution at U+U 20 AGeV collision at b = 4.0 fm. 
In z-plane there is a slight difference between tip-tip collision and sphere-sphere collision. 
On the other hand, in body-body collision the extension of the pressure distribution in 
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Figure 10: The behavior of i>2 as a function of the impact parameter in all 
cases. In this figure elO-b means that the incident energy is 10 AGeV in 
body-body collision and elO-s means that the incident energy is 10 AGeV in 
sphere-sphere collision and e20-t means that the incident energy is 20 AGeV 
in tip-tip collision and so on. 
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x direction is larger than in other cases. Therefore Vi in body-body collision is less than 
one in sphere-sphere and tip-tip collision at any impact parameter. 

Figure 13 shows the pressure distributions at U+U 10 AGeV and U+U 20 AGeV. Since 
there is a slight difference between U+U 10 AGeV and U+U 20 AGeV in the z-plane, the 
growth of V2 in U+U 10 AGeV collision is similar to U+U 20 AGeV. In the y-plane the 
effect of spectator in U+U 20 AGeV is smaller than U+U 10 AGeV. Because the effect 
of Lorentz contraction becomes large in the large incident energy, the spectator becomes 
thin and the influence of it becomes small. Consequently the peak of v% moves to the 
large impact parameter at large incident energy. 

Here, the value of V2 which is influenced by deformation is about 0.015 from fig. |l0[ 
According to the analyses of the excitation function of the directed flow in ref. [21, p2j ], 
the flow becomes slow under phase transition, because the speed of sound becomes small 
in mixed phase, but the effect of the phase transition on flow is expected to be small in 
the AGS energy region. Therefore the effect of deformation is significant and the value of 
V2 which is obtained in this section is important. Further detailed analysis on the relation 
between the effect of phase transition and deformation on flow may provide a new possible 
experimental probe of the phase transition. 
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4 Summary 



We present (3+l)-dimensional relativistic hydro dynamical model of the Lagrangian hy- 
drodynamics without assuming symmetrical conditions. Our algorithm is based on the 
entropy conservation law and the baryon number conservation law explicitly. In our al- 
gorithm we trace the volume elements of fluid along the stream of flux. By using our 
relativistic hydro dynamical model based on the Lagrangian hydrodynamics, the path of 
the each volume element in the phase diagram is able to be traced quite easily. There- 
fore we can investigate directly how the phase transition takes place and affects physical 
phenomenon in an ultra-relativistic nuclear collision. 

Using this model, we have investigated the effect of anisotropic flow in deformed ura- 
nium collisions. The behavior of the flow depends on the pressure distribution and shad- 
owing. Especially shadowing effect increases with the impact parameter. As for differ- 
ences in collision types, there exists only a slight difference between tip-tip collision and 
sphere-sphere collision. On the other hand, in body-body collision is different from 
sphere-sphere collision, i.e. absolute value of t>2 is maximum at b = fm. In light of the 
effect of the incident energy on 112, the peak of V2 is shifted to large impact parameter 
region, because the shadowing effect decreases with increasing the incident energy. We 
accordingly conclude that body-body collision is promising for studying the nuclear equa- 
tion of states and the property of QCD phase transition. These results are consistent with 
ref. H|]. 

We shall remark the following things: In actual experiment, U+U collision is the 
superposition of different types of collision like tip-tip, body-body and so on. Since the 
contribution of body-body collision is not negligible, v% does not equal to zero at b = 
fm. Furthermore, the shadowing has a large effect on v-i- 

In this paper we have applied our relativistic hydrodynamical model to the anisotropic 
flow. Many kinds of applications of our relativistic hydrodynamical model can be con- 
sidered. Using our hydrodynamical model, we can analyze directly the phenomena which 
are sensitive to the phase transition. For example, we can argue how the phase transition 
has effect to the minimum of the excitation function of the directed flow, discussing the 
trajectory of volume element of fluid in T-[i plane. Whether "nutcracker" phenomena 
can be observed in our Lagrangian hydrodynamical model is also interesting problem. 
Hanbury Brown-Twiss effect (HBT) and the influence of anisotropic flow to HBT can also 
be investigated. 

Acknowledgment The authors would like to thank prof. Osamu Miyamura for his 
encouragement to pursue this project and fruitful discussions. They thank Dr. Nobuo 
Sasaki for providing URASiMA and useful comments. They are grateful to prof. At- 
sushi Nakamura for reading the manuscript carefully. Calculations have been done at the 
HSP system of the Institute for Nonlinear Sciences and Applied Mathematics, Hiroshima 
University. This work is supported by the Grant in Aide for Scientific Research of the 
Ministry of Education and Culture in Japan. (No. 11440080) 



References 



19 



For example, see the proceedings of Quark Matter '99, Nucl. Phys. A661, 3c (1999) 

J. Barrette et al. : E877, Phys. Rev. C 56, 3254 (1997) 

C. Pinkenburg et al. : E895, Phys. Rev. Lett. 83, 1295 (1999) 

H. Appelshauser et al. : NA49, Phys. Rev. Lett. 80, 4136 (1998) 

A. M. Poskanzer and S. A. Voloshin for NA49 Collaboration, Nucl. Phys. A661, 341c 
(1999) 

M. M. Aggarwal et al. : WA98, |nucl-ex/980700l 



D. H. Rischke et al. , Nucl. Phys. A595, 346 (1995) 

P. Danielewicz, Roy A. Lacey, P.-B. Gossiaux, C. Pinkenburg, P. Chung, 
J. M. Alexander, and R. L. McGrath, Phys. Rev. Lett. 81, 2438 (1998) 

H. Sorge, Phys. Rev. Lett. 82, 2048 (1999) 

J. D. Bjorken, Phys. Rev. D 27, 140 (1983) 

Y. Akase, M. Mizutani, S. Muroya, and M. Yasuda, Prog. Theor. Phys. 85, 305 
(1991); S. Muroya, H. Nakamura, and M. Namiki, Prog. Theor. Phys. Suppl. No. 120, 
209 (1995) 

T. Ishii and S. Muroya, Phys. Rev. D 46, 5156 (1992) 

J. Alam, S. Raha, and B. Sinha, Phys. Rep. 273, 243 (1996); J. Alam, D. K. Srivas- 
tava, B. Sinha, and D. N. Basu, Phys. Rev. D 48, 1117 (1993) 

U. Ornik, M. Pliimer, B. R. Schlei, D. Strottman, and R. M. Weiner, Phys. Rev. C 
54, 1381 (1996) 

J. Sollfrank, P. Huovinen, M. Kataja, P. V. Ruuskanen, M. Prakash, and R. Venu- 
gopalan, Phys. Rev. C 55, 392 (1997) 

C. M. Hung and E. Shuryak, Phys. Rev. C 57, 1891 (1998) 
B. R. Schlei and D. Strottman, Phys. Rev. C 59, R9 (1999) 

K. Morita, S. Muroya, H. Nakamura, and C. Nonaka, Phys. Rev. C 61, 034904 (2000) 

D. Teaney and E. V. Shuryak, Phys. Rev. Lett. 83, 4951 (1999) 

P. F. Kolb, J. Sollfrank, and U. Heinz, Phys. Lett. B 459, 667 (1999); P. K. Kolb, 



J. Sollfrank, and U. Heinz, [hep-ph/000612S 



D. H. Rischke, Y. Pursiin, J. A. Maruhn, H. Stocker, and W. Greiner, Heavy Ion 
Phys. 1, 309 (1995) 

J. Brachmann, S. Soff, A. Dumitru, H. Stocker, J. A. Maruhn, and W. Greiner, Phys. 
Rev. C 61, 024909 (2000); J. Brachmann, A. Dumitru, H. Stocker, and W. Greiner, 
|nucl-th/9912014l 

E. V. Shuryak, Phys. Rev. C 61, 034905 (2000) 
Bao-An Li, Phy. Rev. C 61, 02193(R) (2000) 
Ben-Hao Sa and An Tai, |nucl-th/991202^ 



20 



[26] C. Nonaka, N. Sasaki, S. Muroya, and O. Miyamura, Nucl. Phys. A661, 353c (1999) 

[27] D. H. Rischke, M. I. Gorenstein, H. Stocker, and W. Greiner, Z. Phys. C 51 (1991) 
485 

[28] Particle Data Group, Eur. Phys. J. C 3, 1 (1998) 

[29] L. D. Landau and E. M. Lifsitz, Fluid Mechanics, Pergamon Press., Oxford(1989) 

[30] M. Gyulassy, K. Kajantie, H. Kurki-Suonio, and L. Mclerran, Nucl. Phys. B237, 477 
(1984) 

[31] S. Muroya and C. Nonaka, Bull. Tokuyama Women's College 6, 45 (1999); |nucl- 
th/9709004| 

[32] A. Bohr and B. Mottelson, Nuclear Structure, Vol. II, P. 133 (Benjamin, New York, 
1975) 

[33] S. Date, K. Kumagai, O. Miyamura, H. Sumiyoshi, and X. Z. Zhang, J. Phys. Soc. 
Japan. 64, 766 (1995) 

[34] N. Sasaki and O. Miyamura, Prog. Theor. Phys. Suppl. 129, 39 (1997) 

[35] N. Sasaki, in preparation 

[36] F. Cooper and G. Frye, Phys. Rev. D 10, 186 (1974) 
[37] U. Heinz, Nucl. Phys. A638, 367c (1998) 

[38] Jean-Paul Blaizot and Jean- Yves Ollitraut, in Quark- Gluon Plasma, R. C. Hwa 
ed. (World Scientific, Singapole, 1990) 

[39] J. Ollitrault, Phys. Rev. D 46, 229 (1992); J. Ollitrault, Phys. Rev. D 48 1132 (1993) 



21 



b = Ofm 




-12 -8 -4 4 8 12 
x [fm] 



6 = 3 fm 




-12 -8 -4 4 8 12 
x [fm] 



6 = 8fm 




-12 -8 




-12 -8 -4 4 8 12 
x [fm] 




-12 -8 -4 4 8 12 
x [fm] 




-12 -8 



Figure 11: The impact parameter dependence of the pressure distribution in U+U 10 AGeV 
sphere-sphere collision at t = 20 fm/c. The left part of figures displays the results at z ~ 
fm and the right part of figures displays those at y ~ fm. The highest values of the pressure 
distribution contour are 0.055 GeV/fm 3 for b = fm, 0.05 GeV/fm 3 for b = 3 and 8 fm, 
respectively and the contour lines are drawn in steps of Ap = 0.005 GeV/fm 3 . 
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Figure 12: The pressure distribution in tip-tip, sphere-sphere and body-body collision U+U 
20 AGeV at t = 20 fm/c. The impact parameter is 4 fm in all figures. The highest value of 
pressure distribution contour is 0.08 GeV/fm 3 for all figures and the contour lines are drawn 
in steps of Ap = 0.01 GeV/fm 3 . 
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E = 10 AGeV 
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Figure 13: The incident energy dependence of pressure distribution in U+U 10.0 AGeV and 
U+U 20.0 AGeV at t = 20 fm/c. The highest value of the pressure distribution is 0.05 GeV/fm 3 
in 10.0 AGeV collision and 0.07 GeV/fm 3 in 20.0 AGeV collision. The contour lines are drawn 
in steps Ap = 0.005 GeV/fm 3 in 10.0 AGeV collision and in steps Ap = 0.01 GeV/fm 3 in 20.0 
AGeV collision. 
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